Anomalous energy transport in the FPU-/? chain 



We consider the energy current correlation function for the FPU-/? lattice. 
For small non-linearity one can rely on kinetic theory. The issue reduces 
then to a spectral analysis of the linearized collision operator. We prove 
thereby that, on the basis of kinetic theory, the energy current correlations 
decay in time as t~ 3 / 5 . It follows that the thermal conductivity is anomalous, 
increasing as N 2 / 5 with the system size N. 

1 Introduction and physical background 

With the availability of the first electronic computing machines, Fermi, Pasta, and 
Ulam [1 j investigated the dynamics of a chain of nonlinear oscillators, in particular, 
their relaxation to thermal equilibrium. Their work had a, in retrospect surprisingly, 
strong impact. We refer to the special issue [2] which accounts for the first fifty 
years. In our contribution, we will study the /3-chain. This is a linear chain of 
equal mass particles which are coupled to their nearest neighbors by nonlinear 
springs with a potential of the form Up(r) = |r 2 + |/3r 4 , with (5 > 0, and r 
being the string elongation. (According to the FPU convention, the a-chain has a 
nonlinearity ^ar 3 instead of \(3r A .) If we denote the particles positions by G M, 
and their momenta (velocities) by p, £ M, then the /3-chain has the Hamiltonian 
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Pi = U'p(q i+ i - qi) - Upfai - qi-i). 
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Over the last decade there has been a lot of interest to understand the energy 
transport through one-dimensional chains, amongst them the FPU /3-chain [3 ]. Nu- 
merically, one common setup is to consider a chain of length N, and to couple 
its left- and rightmost particles to thermal reservoirs at temperatures T_ and T+, 
respectively. For long times the chain relaxes to a steady state with a non-zero 
average energy current j c {N) = (T_ — T + )./V _1 k(./V), and the interest lies in the 
dependence of k(N) on N for large N. For a regular transport, i.e., for trans- 
port satisfying Fourier's law, one has k(N) — ► const, for large N. Anomalous 
transport corresponds to k(N) ~ N a , with < a < 1. In the /3-chain more 
recent molecular dynamics simulations point to an a of approximately 0.4 El HI, 
and thus a larger energy transport than expected on the basis of Fourier's law. In 
these simulations chain lengths of up to N = 2 16 are used, and the result seems 
to be stable for a range of fairly low boundary temperatures. In [61 it is claimed 
that for somewhat higher boundary temperatures, there is a crossover at large N to 
k(N) ~ N3. Hence, even on the numerical level the accurate value of a is still 
being debated. 

In this paper, we will adopt a different, but physically equivalent procedure. 
One prepares initially (for t = 0) the infinite /3-chain in thermal equilibrium at 
temperature T > 0. This means that the initial conditions of the Hamiltonian 
dynamics are distributed according to the (at this stage formal) Gibbs measure 

Z-^- H/T II [d*d»] • (1-3) 

This measure does not change in time. One now adds some extra energy close to 
the origin and studies the spreading of this excess energy. To be more precise, let 
us introduce the local energy, e^, at the site iGZby 

ei(q,p) = - [pj + Ufs{q i+ i - + Up{qi - . (1.4) 

We also employ the shorthand notation ei{t) = ei(q(t),p(t)), where (q(t),p(t)) is 
the solution to the Hamiltonian dynamics (11.21 ) for given initial conditions. Then 
we define the normalized local average excess energy by 

S(i,t) = - ((ei(t)eo(O)) - <ei)(eo)) . (1.5) 
X 

Here (•) denotes the thermal average (11.31 ) over the initial conditions, and \ = 
J2i (( e i e o) — ( e i)( e o)) is a normalization guaranteeing £)ii?(*>0) = 1. One has 
S(—i, t) = S(i, t), and the energy spread at time t is defined as the spatial variance 

D(t) = ^2i 2 S(i,t). (1.6) 

Fourier's law corresponds to a diffusive spreading, D(t) = 0(t) for large t, 
while an exponent a > corresponds to superdiffusive spreading with D(t) = 
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0(t 1+a ). These properties can be more conveniently reformulated by introducing 
for each directed bond from i to i + 1 a current ji^+i, so that the energy continuity 
equation holds in the following form: 

+ - ji-ij = 0. (1.7) 
For the FPU-/3 model such a current observable is given by 

ji,i+i{q,p) = -^{Pi+i +Pi)U'p(q i+ i - qi). (1.8) 

Obviously, (ji^+i) = 0. We next introduce the energy current-current correlation 
function 

C/#) = £<io,i(t)ji,m(0)>. (1-9) 

Then 

D{t) = D(0) + - f ds I ds'Cp(s-s'). (1.10) 
X Jo Jo 

Note that \Cp(t)\ < Cp(0) < oo. Clearly, if / °° dt\Cp(t)\ < oo, then D(t) = 
0{t). On the other hand, if Cp(t) = 0{t a - v ) for large t with < a < 1, then 
D(t) = 0(t 1+a ), and the spreading is superdiffusive. 

The problem of regular versus anomalous energy transport may thus be re- 
phrased as whether Cp(t) decays integrably or not. Unfortunately, such a reformu- 
lation is of little help. To estimate the decay of a time correlation in equilibrium, 
such as Cp(t), is an exceedingly difficult problem. However, in the limit of small 
j3, through methods from kinetic theory, Cp(t) may be expressed in a more acces- 
sible form. For the complete argument we refer to El [8] [9]]. Here we only state the 
small (3 form of Cp(t). To do so will require some preparation. But the goal of our 
contribution is to estimate the decay of Cp(t) for the FPU-/? chain in the limit of 
small [3. 

At P = 0, the system reduces to the harmonic Hamiltonian 

H(q,p) = M + g(«+l " ft) 2 ] , (LID 
which has the dispersion relation 

Lo(k) = y^(l — cos k) 



k 

sin — 
2 



(1.12) 



Here we use the convention that the discrete Fourier transform yields 27r-periodic 
functions, and also declare that the term "periodic function" always refers to a 
function which is 2-7r-periodic in all of its arguments. It will be convenient to choose 
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as the basic periodic cell the interval I = [0, 2tt). In particular, then x mod 2ir G / 
for all On / the dispersion relation is simply 



and thus also for all x ^ 0, 



x 

u(x) = sin-, (1.13) 



1 x 

u (x) = - cos -, (1.14) 



and we let arbitrarily u>'(0) = 0. We also introduce 

Q(x, y, z) = uj(x) + oo{y) - uj{z) - uj(x + y - z) (1-15) 

for x, y, z £ R. With these conventions the linearized collision operator of the 
FPU-/? lattice in the kinetic limit is given by 

(Lf)(x) = Jdy Jdz5(n(x,y,z))(f(x)+f(y)-f(z)-f(x + y-z)), 

(1.16) 

with / periodically extended from / to R, see [8]. 

L describes the collision of two phonons, where x,y label the incoming mo- 
menta and z, x + y — z label the outgoing momenta, thus by fiat satisfying mo- 
mentum conservation modulo 2n. Through the 5-function the collisions are also 
constrained to conserve energy. Note that at this stage, the definition in (11.161 ) is 
only formal since no prescription is given of how to deal with the 5-function. It 
turns out to be useful to consider L = luLuj as a linear operator on L 2 (I), with oj 
being understood as the multiplication operator by the function oj. We will prove 
later that L is a bounded positive operator with a decomposition 

L = W-A, (1.17) 

where A is compact and W is a multiplication operator. 

Now we are in a position to state the conjectured behavior of Cp(t) for small 
coupling f3. 



cA (1.18) 



Kinetic conjecture: For any iGR and temperature T > 

T 2 i r 
lim Cfl(/3~ 2 t) = — (w',exp -7r _1 (12T) 2 |t|£ 

/3^0+ Z7T \ L 

where (-, ■) denotes the scalar product in L 2 (I). (A more detailed discussion about 
the scaling factors can be found in [8].) 

Thus for small (3, the decay of Cp(t) is obtained from the spectral properties of 
L, certainly a more accessible item than the full Hamiltonian dynamics. Our goal 
here is to study the behavior of the kinetic correlation function 

C(t) = <a/,e-l*IV). (1.19) 
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In kinetic theory, it is a common practice to use the relaxation time approxima- 
tion, which in our case amounts to dropping the operator A, that is, to approximate 

(u/,e-l'IV) « (u',e-^ w oj'}. (1.20) 

As we will show, W(x) = W(2ir — x), and for < x <C 1, W(x) behaves asymp- 
totically as x 5//3 . Thus the relaxation time approximation predicts (a/, e~^ L u)') = 
0(t~ 3 / 5 ) for large t, as has been derived in ifTOl . 

L has the range of W as its essential spectrum. In particular, the essential 
spectrum starts from 0. Thus it is not obvious that the asymptotics predicted by 
the relaxation time approximation is really the correct one. To understand the time 
decay leads to two distinct mathematical issues. 

(1) The so called collisional invariants, which in essence are zero modes of L, 
could in principle prevent C (t) from decaying to 0. To exclude such a possibility, 
we have to characterize all collisional invariants, which involves solving a non- 
trivial functional equation. 

(2) We will use the resolvent expansion to estimate (a/, e~^ L uj'). In our case, it 
turns out to be necessary to make the expansion to the second order, yielding 



A + L / \ X + W I \ ' X + W X + W 

+ T77 A C 1 - 21 ) 

\ X + W A + L X + W I 

The first term is identical to the relaxation time approximation, and behaves as 
A~ 2 / 5 for < A < 1. The second and third term will be shown to be C(A _1 / 5_£ ) 
for any e > 0. Although also this second contribution is divergent, the first term is 
dominant, and thus we confirm the prediction of the relaxation time approximation 
in this particular case. 

An inherent difficulty in resolvent expansions is the estimation of the remainder 
term, such as the last term in (I1.21I ). Our method bears some similarity to the 
Birman-Schwinger estimates used in quantum mechanics. It relies on the fact that 
the resolvent expansion is made up to an even order, as well as on the operator 
B = W~ l l 2 ' AW~ X I 2 being compact. In fact, it is likely that similar techniques 
can be used to study many of the cases where a decomposition L = W — A, with 
W > and a compact B, is possible, although we would expect the optimal order 
for the resolvent expansion to vary from case to case. The exact order, as well as 
the exact power of the decay, would naturally depend also on the function a/. A 
reader interested in such generalizations is invited to jump ahead to the proof of 
the main theorem in Section [6] 

Our results imply that, on the kinetic time scale, the energy spread is superdif- 
fusive, with D(t) ~ ct 7 ^, c > 0, for large t. This corresponds to a heat conduc- 
tion exponent a = | and is in agreement with the molecular dynamics simulations 
of 6J0. As the example of long time tails in classical fluids teaches us, kinetic 
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theory might miss the true asymptotic decay of equilibrium correlation functions. 
Whether this is the case also for the FPU-/5 chain, remains a challenge for the 
future. 

From the point of view of kinetic theory, our result is fairly surprising. Usually 
linearized collision operators have a spectral gap implying exponential decay of the 
current-current correlation function, and diffusive spreading for the corresponding 
conserved quantity. In fact, we are not aware of any other Boltzmann type kinetic 
model which would exhibit superdiffusive spreading. 
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2 Main results 

To define L, we first need to find all solutions to the energy constraint. The solution 
manifold to U(x, y, z) = 0, is clearly non-empty, as there are the trivial solutions 

z = x and z = y. (2.1) 

We will later prove in Corollary I3.3l that. in addition, there is a solution y = h(x, z), 
and that all other solutions are modulo 2ir equal to one of these three. For x,z € I 
the function h is given by 

, , . z — x . ( \z — x\ x + z\ 

nix, z) = h 2 arcsm tan cos (2.2) 

V ; 2 V 4 4 / 

where arcsin denotes the principal branch with values in [— n/2, vr/2]. We extend 
h to R 2 by defining 

h(x, z) = h(x mod 2tt, z mod 27r) — i(x), (2.3) 

where i(x) = x — (x mod 2tt) £ 27rZ. This choice makes h everywhere continu- 
ous while ensuring that for all x, z G R, we still have h(x, z), z) = 0. 

The energy conservation ^-function can then be formally resolved by integrat- 
ing over some chosen direction: for instance, choosing the y-integral for this pur- 
pose would yield for any z/x and for any continuous periodic function G, 

dy d(Q{x,y, z))G(x,y, z) 
1 



G(x,z,z) + ——— r — rvG(x, h(x, z),z). (2.4) 



\d2&(x, z, z)\ ' ' \&2^l(x,h(x,z), z)\ 
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However, this procedure is somewhat suspect here, as it will lead to terms of the 
type oo — oo, related to the trivial solutions and canceled only due to symmetry 
properties. An additional difficulty lies in the application of the definition to func- 
tions G which are not continuous but merely L 2 -integrable. To put the definition 
of L on a firmer ground, we will resort to a different approach in Section [3j we 
replace 5 in (11.161 ) by a regularized (5-function S e (X) = e7r _1 (e 2 + X 2 )~ l , e > 0, 
and then show that there is a unique self-adjoint operator L which agrees with these 
operators in the limit e — > 0. Our choice of regularization for the (^-function is not 
completely arbitrary: in the kinetic limit of lattice systems with random mass per- 
turbations the corresponding 5-f unction also appears via a sequence of S e -functions 
(see, for instance, Proposition A.l in ifTTTO . 

A somewhat lengthy computation, to be discussed in Sections [3] and 01 shows 
that the formal procedure explained before is essentially correct: the trivial solu- 
tions give no contribution, and the unique limit operator L is 

L = V + K 1 -2K 2 , (2.5) 

where K\ and K 2 are integral operators determined by the integral kernels 

K 1 (x,y) = A 1{F -l X >y )> ° ) and K 2 (x,y) = J (2.6) 

which are defined for x,y E I using the auxiliary functions 

x y\ 2 ... x . y 



F±{x,y) = [cos- + cos-j ±4sin-sin-. (2.7) 
In addition, V denotes a multiplication operator by the function 

V(x) = jdyK 2 {x,y). (2.8) 

L was already used as the linearized collision operator in [10]. In addition to L, 
L = ujLoj, and W = u 2 V, the operator B = W- l / 2 {W - ^W' 1 / 2 will be of 
importance. Explicitly, B is then defined via the integral kernel 

B(x,y) = V(xr 1/2 (2K 2 (x,y) - K x (x, yWiy)" 1 / 2 . (2.9) 

Let us next list the main properties of these operators, to be proven in Sections 
0]and[5] We start with the results related to item Q] mentioned in the introduction. 

Definition 2.1 A measurable periodic function ip : R — > C is called a collisional 
invariant if for almost every x,y,z G R such that £l(x,y,z) = 0, 

rp(x) + rp(y) - ip(z) -ip(x + y- z) = 0. (2.10) 

In the definition, "almost every" refers to the Lebesgue measure on any two- 
dimensional submanifold of the full solution set. The following theorem shows 
that, in the case considered here, there are only the obvious collisional invariants. 
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Theorem 2.2 Suppose tp is periodic and locally integrable: G L l (I). Then 
ij) is a collisional invariant if and only if there are C\,C2 £ C such that ip(x) = 
c\ + C2Uj{x) for a.e. x. 

In higher dimensions there is a general argument which identifies the collisional 
invariants under minimal assumptions on oj lfl2l . In contrast, our proof here relies 
heavily on the specific form of u, and does not exclude the appearance of non- 
trivial collisional invariants in some other one-dimensional systems. 

Definition 2.3 We define a parity transformation P : L 2 (I) — > I? (I) by letting 
(Pi/;)(P) = i/>(0) and, for x 6 (0, 2vr), 

(Pip)(x) =ip(2ir - x). (2.11) 

Clearly, u>(x) is symmetric, and lo'{x) is antisymmetric under P. 

Proposition 2.4 L is a bounded positive operator, and B is a compact self-adjoint 
operator on L 2 (I). Both B and L commute with the parity operator P. In addition, 
B < 1, and Bip = ip if and only if the periodic extension o/ V _1 / 2 ?/> is a collisional 
invariant. 

Theorem 2.5 Let R: (0, oo) — ► R + be defined by 

tf(A) = (t/,^-j-W). (2.12) 

Then there is < cq < oo such that with a = | 

lim \ a R(\) = co. (2.13) 

Since R(X) = / °° dt e~ xt C(t), for A > 0, R(X) is a Laplace transform of the 
monotonically decreasing positive function C(t). Methods from Tauberian theory 
can then be used to connect the asymptotic behavior of R and C, proving that the 

3 

asymptotic decay of the current-current correlations is given by C(t) = 0(t~5 ), 

t 2 

and that the integrated correlations grow like y Q dsC(s) = 0{U). Explicitly, 

Corollary 2.6 With cq > and a = | as in Theorem Xl. 51 and with V denoting the 
gamma function, we have 

lim t^Cit) = (2.14) 
t^oo 1 (a) 

(For a proof of the result, see for instance "Zusatz zu Satz 2" on p. 208 of lfT3l .) 

We have divided the proof of the above results in four sections. We solve 
the energy constraint and derive the above form for the operator L in Section [3] 
Proposition I2.4l is proven in Section |4l which includes, in particular, the estimates 
proving the compactness of B. We study the collisional invariants in Section \5\ 
and prove Theorem 12. 21 there. Finally, these results are then applied in a resolvent 
expansion, and we prove Theorem 12.51 in Section [6l The short Appendix contains 
a convenient estimate for the norm of an integral operator. 
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3 Resolution of the energy constraint 



We will define the operator L by the following procedure: we consider a regular- 
ization of the 5-function by 

8e(X) = -^—2, (3.1) 

for any < e < 1. Let L e denote the operator defined by (11.161) after 5 has been 
replaced by 5 e . This yields a bounded operator, for which using the symmetry 
properties of the integrand 

i(f,L e f)= [ dxdydz5 € (n(x,y,z))\f(x) + f(y)-f(z)-f(x + y-z)\ 2 . 

(3.2) 

Our aim in this section is to prove the following result about the limiting behavior 
of this quadratic form when e — > + . 

Proposition 3.1 For any f : R — > C, which is periodic and Lipschitz continuous, 
the limit lim e _ )0 + (/; L e f) exists, and it is non-negative, finite, and equal to 

dxdz 1 |/(x) + f(h(x, z)) - f(z) -f(x-z + h(x, z))\ 2 
P 2 y / F + (x,z) 

dxdzf{x)*{V{x)+K 1 (x,z)-2K 2 {x,z))f{z). (3.3) 

The proof of the Proposition will require some fairly technical estimates not needed 
later, and a reader accepting our definition of the operator L and the equality in 
(13.3b can well skip the proofs of the Lemmas below in the first reading. 

We will begin by constructing the solutions to the energy constraint O = 0, and 
then study the behavior of Q around this set, in order to evaluate the limit of the 
approximate ^-functions. Let D = [0, 2n] 3 be the closure of I 3 . For (x, y, z) £ D, 



tt(x, y, z) = sm - + sm - - sm - 



. x + y - z 
sm 



(3.4) 



Since then — tt < a+ | z < 2tt, we can split D into two sets U+ and £/_, depending 
on the sign of the last term. Explicitly, 

U + = {(x, y,z)eD\x + y-2ir<z<x + y}, (3.5) 
f/_ = {(x,y,z) G D | x + y < z or z < x + y — 2tt} , (3.6) 

and fl(x, y, z) = fi, a (x, y, z) with a = +1 if (x, y, z) G U+, and with a = —1 if 
(x, y, z) G where 

x y z x + y — z 
n a (x,y,z) = sin - +sm- - sin - - asm . (3.7) 
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The following representations of these functions will become useful later (they can 
be checked, for instance, by expressing the trigonometric functions in terms of 
complex exponentials): for all x, y, z G R, 

n/ \ a • x ~ z ■ y~ z ■ x + y n Q^ 

n + (x,y,z) = 4 sin—— sin—— sin—— , (3.8) 

„ , s „ / x + z . x — z x — z . 2y + x — z\ 

il-(x, y, z) = 2 cos sin h cos sin . (3.9) 

V 4 4 4 4 / 

From these, we directly find the zeroes of S7: 

Lemma 3.2 Lef Z = {(x, y, z) G -D | y, z) = 0}, and denote Z± = Z n f/±. 
Z + consists of those (x,y,z) G D for which either z = x, z = y, or y = x G 
{0, 2-7r}. Z_ consists of those (x, y,z) G D which satisfy any of the following three 
conditions, where h is defined by (12. 21) and (12. 3\) . 

1. x = 0, z = 2ir, or x = 2ir, z = 0, 

2. x < z, and y = h(x, z), 

3. x > z, and y = 2tt + h(x, z). 

In addition, for (x, y, z) G J7_, w/?/z x ^ z, we have sign(z — x)d y Q-(x, y, z) > 
cos 2 

Proof: By (13.81 ). ^ = on £/+ if and only if one of the three factors is zero. Since 
l^f^l < f , and < ^ < it, for any a,b G {x,y,z}, this can be checked to 
coincide with the above classification of Z + . 

To compute Z_, let us first consider the case cos = 0. Then either x = 0, 
2 = 27r, or 2 = 0, x = 2tt, and both cases can be checked to form solutions for 
any y. Otherwise, cos ^J 5 > 0. Also cos ^J 2 > cos ^j 2 , as < < tt. 

Similarly, as |£=2| < 27r ~ a: + 27r - 2 < 7r> we have cos ^J 2 > — cos ^x 2 . Therefore, 
| cos < cos Also by (f379b 

cJyO-^, y, z) = cos ^— — cos - ^ — -. (3.10) 

We then split the proof into three steps with additional conditions onx, z. 

Assume first z = x. Then (x, y, z) G U- if and only if y = or y = 2n, and 
both cases clearly yield solutions. Since h(x, x) = 0, both cases are covered by 
the Lemma. 

Assume then z > x. Then (x, y, z) G U- if and only HO < y < z — x < 2ir. 
This implies that 1 2 v + ^~ z | < < |, and thus in this case d y Q- > cos 2 ^f 2 > 
0. On the other hand, by explicit computation, then f2_ (x, 0, z) < and 0_ (x, z — 
x, z) > 0. Therefore, for such x, z there is a unique solution y G [0, z — x]. By 
(13.91) this solution satisfies 

. 2y + x-z x + z z-x 

sin = cos — - — tan — - — . (3.11) 
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This equation has infinitely many solutions !/£R, but the above bounds show that 
exactly one of them, 



z — x 



+ 2 arcsin 



x + z z — x 
cos — - — tan — - — 



h(x,z), 



(3.12) 



with arcsin denoting the principal branch with values in [— 7r/2, 7t/2], can belong 
to [0, z — x\. Since there must be a solution in this interval, we find that h(x, z) G 
[0, z — x], and thus (x, h(x, z),z) G ?7_. 

To complete the analysis, assume z < x. Then (x, y, z) G U- if and only if 
2-k + z — x < y < 2ir. We let x' = 2ir — x, etc., when z' > x', and < y' < z'—x'. 
As always £7_ (x', y', z') = y, z), we can conclude that for any x, z there is 

a unique solution in U- which satisfies y' = h(x', z'), i.e., the solution is 



y = 2vr + 



+ 2 arcsin 



x + z x 
cos — - — tan — 



2it + h(x,z). (3.13) 



It also follows that in this case, d y £l(x, y,z) < — cos < for all 2ir + z — x < 



y < 2ir. This completes the proof of the Lemma. 



□ 



As f2 is periodic, the Lemma yields immediately also a classification of the zeroes 
of tt in M 3 . 



Corollary 3.3 f2(x, y, z) =0 if and only if at least one of the following equalities 
holds modulo 2ir: z = x, z = y, or y = h(x, z). 

Proof: It is clear from the Lemma that any (x,y,z) G M 3 satisfying the above 
condition is a zero of Q. For the converse, assume y, z) = 0. Then for 
x' = x mod 2tt, etc., also y', z') = 0, and we can apply the Lemma. If 
(x 1 , y', z') G U + , then either z' = x' , z' = y' or x' = y' = 0. Since the last con- 
dition implies z' = 0, and thus y 1 = = h(0, 0), also the last instance is covered 
in the Corollary. If the point belongs to U-, we must have y' = h(x', z') mod 27r, 
and thus then y = h(x, z) modulo 27r. □ 

The following Lemma can then be used to compute the relevant limits for inte- 
grals over U- : 

Lemma 3.4 Suppose G : R 3 — > C is a periodic continuous function. Then 
lim / dxdydz 5 e (Q(x, y, z))G(x, y, z) 

= I dxdz - 2 =G(x,h(x,z),z). (3.14) 
Jp v / F + (x,z) 

Proof: As G is periodic and continuous, it is also bounded. Let < e < 1 be 
arbitrary, and denote X £ = [0,e] x {2ir — e,2ir] U [2ir — e,2ir] x [0,e]. Let us 
first consider some fixed x,z G [0, 2ir] 2 \ X £ , x ^ z. By Lemma l3~2l \d y £l\ > 
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cos 2 > sin 2 | > 0, and y h-> Q[x,y,z) is a bijection with a unique zero, 
yo, which is equal to h(x, z) modulo 2ir. We change the integration variable y to 
t = Q,(x, y, z)/e, which shows that the integral over y is equal to 

*>A 1 1 

a/e \d 2 VL{x,y{et),z)\Tt(l + t 2 ) 

with a < and b > and y(0) = yo- This is always bounded by a constant which 
depends on e but not on e. Thus an application of the dominated convergence 
theorem shows that 



lim / dxdydzt((x, z) X e ) 8 e (Q,[x, y, z))G(x, y, z) 

«-o+ Ju- 

dxdz .^ X / z) J X ^ G(x, h(x, z),z). (3.16) 
p \o2\l-[x,h(x,z),z)\ 

We used here the observation that the set of x = z, as well as that of (x, z) for 
which a = or b = 0, have zero measure. Here, by (13.101 ). we have 



\d2&-(x, h(x, z), z)\ 



x — z 

COS ; 



. 2 2h + x - z \ 2 
1 — sin ■ 



1 

2 X Z 2 X Z 2 X ~\~ Z \ 2 1 



cos — sm — ■ — cos 

4 4 4 



V = ~y/F+(x,z), (3.17) 



where the last equality can be checked by a calculation, for instance, using the 
identity cos 2 u = |(1 + cos(2n)). For all x, z G /, we clearly have an estimate 

0<K 2 (x,z)= 2 < (sin^sin^)" 172 . (3.18) 

— 1/2 

Thus F + is integrable, and we can again apply dominated convergence to take 
the limit e —> inside the integral. This proves that the right hand side of (13.161) 
converges to the right hand side of (13- 14b . 

Therefore, to complete the proof of the Lemma, it is sufficient to prove that the 
integral over (x, z) G X £ vanishes when first e — > and then e —> 0. In fact, using 
the symmetry between the two components of U- and the boundedness of G, it is 
sufficient to study the integral 

re r2Tr rz—x j 

J £ = dx dz dy- (3.19) 
Jo h-K-e Jo Tre 2 + n± 

We split the integral over y into two parts at y = 7r. If < y < it, we have 
2d z £l- > —cos | > cos|. Therefore, we can perform the z integral first, as 
above, and conclude that the result of the z-integral is uniformly bounded in e and 
e. Performing then the x and y integrals, shows that the full integral is bounded by 
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a constant times e. In the remaining region, ir < y < z — x, we have 2<9 a .O_ > 
cos | > cos |. Thus in this case, we can perform the x integral first, with a 
uniformly bounded result. Performing then the y and z integrals, and combining 
the bound with the earlier estimate, proves that there is c > such that J £ < ce. 
Thus J £ —> as e — » 0, which concludes the proof of the Lemma. □ 

To complete the proof of Proposition [3j] we require one more Lemma, closely 
related to the above estimates. 

Lemma 3.5 Assume G is measurable and periodic on IR 2 . Then 

I dxdz J_^ G(x,h(x,z))= [ dxdy2 HF -^l^ G(x,y), 

(3.20) 

as long as either G > 0, or one of the above integrals is absolutely convergent. 

Proof: As is apparent from (13.20I ). the proof is accomplished by a change of inte- 
gration variables from z to y = h(x, z) for a fixed x. However, even computing 
the local inverse functions from (12.21 ) does not appear to be completely straightfor- 
ward. We will resort to a roundabout way, which relies on the fact that h(x, z) is a 
solution to the energy constraint on fL_. 

If h(y;x) is a local inverse of h(x,-), then for all y in its domain there is 
n G {0, 1} such that (x, y + 2irn, h(y; x)) G U- and 

tt(x,y,h(y;x)) = n„(x,h(x,h(y;x)),h(y;x)) = 0. (3.21) 

Conversely, assume that x G (0, 2ir) is given, and h(y; x) is a map from some 
interval J C (0, 2ir) to / such that either x + y < h(y; x) < 2ir or < h(y; x) < 
x + y — 2ir for all y, and Q(x, y, h(y; x)) = 0. Then by Lemma [3T2l we must have 
in the first case y = h(x, h(y; x)), and in the second case, y = 2ir + h(x, h(y; x)) 
for all y. 

Therefore, to construct all possible local inverse functions of h(x, •), it is suffi- 
cient to find for given x, y all z such that (x, y, z) G t/_, and f2_ (x, y, z) = 0. We 
begin from the following representation of for all x, y, z, 

( x + y x — y x + y x + y — 2z\ 
\l-(x,y,z) = 2 sm — - — cos — h cos — - — sm . (3.22) 



Let us assume that (x,y, z) G U- with < x < 2ir. Then cos = implies 



4 

y = 2tt — x, and thus then f2_ = 2 sin § > 0. Thus if Jl_ = 0, we have 



. x + y — 2z x + y x — y 
sm = — tan cos . (3.23) 
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Since x, y, z are real, this is possible only if the absolute value of the right hand 
side is less than or equal to one. This condition is equivalent to the condition 
F(x, y) > 0, with 

2 x + y 2 x + v 2 x ~y 

F(x,y) = cos — sm — - — cos — - — . (3.24) 

A brief computation reveals that, in fact, F(x, y) = \F-.{x, y), and thus F_ > 
is a necessary condition to have any solutions. 

When F_ (x,y) > 0, (13.231 ) holds if and only if there is n G Z such that z = z n , 
where 

O , X + y i t lV»o ■ f + X + y X ~V\ ,~ OCk 

z n = 2im H h (—1) zarcsm I tan — - — cos — - — , (3.25) 



with arcsin denoting the principal branch. There can be maximally two values of n 
for which z n belongs to [0, 2n). However, by inspecting the sign of 8 3 Q-, similarly 
to what was done in the proof of Lemma 13.21 we find that for F_ > there are 
exactly two solutions in U- , and that for a given x there are maximally two values 
of y for which F_ = 0. If F_ > 0, the solutions are explicitly z = h±{y;x), 
where for either choice of the sign a G {±1} 

i i \ x + y i o i o ■ (*. x + y x ~y \ 

n a [y; x) = — h c a 2ir + a 2 arcsin tan — - — cos — - — , (3.26) 



and c a = 0, if a = +1, and c a = {-\)^+v>^), if a = -1. 

Therefore, apart from a finite number of values y G [0, 2ir), there are either no, 
or there are exactly two, solutions in U- . Both of the solutions satisfy (13.231 ) and 
thus also for any such z 



\d 3 n(x,y,z)\ 



x + y 
cos — ■ — 



! ■ 2 x + y - 2z 

1 — sm ■ 



= y/F(x,y) = l - y/F-(x,y). (3.27) 
On the other hand, by implicit differentiation we find 

82^1 ~(x, h(x, z), z)d z h(x, z) + 83ft -(x, h(x, z), z) =0. (3.28) 
By (13.171) this implies 

±y/F+(x,z)\d g h(x,z)\ = \8 3 n4x,h(x,z),z)\ , (3.29) 

which allows to compute the Jacobian of the change of variables. 

Collecting all of the above results together, and applying Fubini's theorem, we 
can conclude that (13.201 ) holds for G > and for any G which is bounded. This 
implies, in particular, that the integrals are equal if G is replaced by \G\ for any 
measurable G. Thus if either of these integrals in (13.201) is absolutely convergent, 
then the other must be so as well. Then an application of dominated convergence 
theorem proves that (13.201 ) holds also for such measurable G. □ 
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Proof of Proposition \3.1[ Let / be a periodic Lipschitz function. We express 
(f,L e f) as an integral over I 3 using (I3.2I ). As D = U + U U-, and D \ I 3 and 
£7+ D t/_ have measure zero, we can split the integral into two parts by using 
J I3 = J v + . Since the factor multiplying 5 e in the integrand is positive, pe- 
riodic, and continuous, Lemma [3~4l implies that the integral over {/_ converges to 
the left hand side of (13.3I ). By boundedness of / and applying (13.181) . the integral 
yields a finite, non-negative result as claimed in the Proposition. 

Thus in order to prove convergence to the left hand side of (13.31 ). we only need 
to show that the integral over U + vanishes as e ^ + . For any (x, y, z) G U + , 
using (13.81) and the fact that | sinx| > ||x| for |x| < §, 

m / m ^ A \ x ~ z \ \y~ z \ ■ x + y ^ 1 , M , / x 

P+{x,y,z)\ > 4—^- ^— sin — > -^\x - z\\y - z\m(x,y), 

(3.30) 

where m(x, y) = min(x + y, 4-7T — x — y). This implies that, if \x — z\ < \y — z\, 
then by the Lipschitz property of /, there is a constant c > such that 

irrWl/OzO + f(y) - /(*) - /(^ + y - ^)l 2 < , c , 2 , ^if u - 

e 1 + m(x, y) 1 (vre) z + \x — z| 4 

(3.31) 

H" \y — z\ < |x — z|, the same estimate holds after x and y have been interchanged 
on the right hand side. Therefore, 

[ dxdydz -i ^|/(x) + /(y) - f(z) -f(x + y- z)\ 2 

Ju+ tt e 2 + il(x,y,zy 

f , , l(m(x,y)>e) f 00 , e t 2 ,212 
<c / dxdy V y tf; ~ 7 / dt- 7 . 2 ■ . +c / e 2 6 1 - 2 , (3.32) 
Jp m(x,y) z 7r(7re)^+i 4 

where the second term is an estimate for the integral over (x, y) with m(x,y) < e 
— these are contained in the two boxes [0,e] 2 and [2ir — e,2ir] 2 and we have 
there estimated the integrand trivially using Q, 2 > 0. The first integral over x, y is 
0(| lne|), and the second integral is 0(e 1 ^ 2 ), as seen by changing the integration 
variable to s = e~ l l 2 t. Therefore, we can conclude that the left hand side of (13.321) 
vanishes as e — ► + , i.e., that the integral over U + does not contribute to limit, as 
long as / is a Lipschitz function. 

Thus to complete the proof of the Proposition, we only need to prove the equal- 
ity in (13.31) . Instead of doing this directly, let us come back to the integral over t7_, 
which was proven above to converge to the left hand side of (I3.3I ). The set U- is 
clearly invariant under x <-> y. By inspection we check that this is also true for the 
map z 1— ► z', with z' = x + y — z + 2tt, for x + y < z, and z' = x + y — z — 2ir, oth- 
erwise. Similarly, f/_ is left invariant under the map (x, y, z) 1— ► (a/, y', z'), with 
y' = z, z' = y, and x' = x + y — z + 2tt, for x + y < z, and x' = x + y — z — 2tt, 
otherwise. All of these maps leave also $7 invariant. We can thus first expand the 
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square and then use the above mappings to appropriately change variables to prove 
that 

~ J dxdydz5 e (n(x, y, z))\f{x) + f(y) - f(z) -f(x + y- z)\ 2 

= [ dxdydz5 e (n(x,y,z))f(x)*(f(x) + f(y) -2f(z)). (3.33) 
JU- 

Lemma 1341 can be applied to the right hand side proving that it converges to 

/ dxdz 2 f(x)*(f(x) + f(h(x,z))-2f(z)) (3.34) 

JP y/F + (x,Z) 

By Fubini's theorem, the first and the last terms are equal to those of the right 
hand side of (13.3I ). The middle term is absolutely convergent by (13- 18b . and thus 
an application of Lemma l3~5l shows that it is equal to the missing K\-te.rm in (13.31) . 
This completes the proof of the Proposition. □ 

4 Linearized collision operator 

We will derive in this section the regularity properties of L and B and prove Propo- 
sition l2.4l Let us recall the definition of L e in Section[3l and Proposition 13 . 1 I pro ven 
there. Let also L e = ujL e uj. If / is a periodic Lipschitz function, then so is g = uf, 
and thus an immediate consequence of the Proposition is that then 

lim -(f,L e f) = lim (g,L e g) = (f,Lf) > 0. (4.1) 

Here we have employed the definition of L to identify it in the right hand side of 
(I3.3I ). We shall soon prove that L is a bounded operator on L 2 (I). As Lipschitz 
functions are dense in L 2 (I), the above result implies that L is a positive operator. 
Moreover, L is then uniquely determined by (14.11 ) in the following sense: Suppose 
L' is another self-adjoint operator (not necessarily bounded) for which (/, L'f) = 
(/, Lf) for every Lipschitz function /. Since then V — L is densely defined with 
(/, (1/ — L)f) = 0, we can conclude using the polarization identity that V f = Lf 
for all / Lipschitz. As L is bounded and self-adjoint, this implies L' = L. Thus 
we only need to check that L is bounded, and to show that [P, L] = 0, in order to 
conclude the properties stated about L in Proposition 12.41 We remark in passing 
that only L will be proven to be bounded, the operator L could well be unbounded. 

For the proof of compactness of B, we need more precise estimates on W 
and on the kernel functions K\ and K2. We recall the definition of W given in 
SectionEl W = uj 2 V. As to be shown, the exponent a = | in the main theorem is 
determined by the behavior of W(x) near x = 0. This will be summarized in the 
following Lemma. 
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Lemma 4.1 The function W : K — » w symmetric, PW = W, and continuous. 
In addition, there are constants C\,C2 > 0, such that for all x GM 



ci 



x 

8111 — 

2 



x 

sin — 
2 



aw(i also lim x ^o ( | sin ■ 



< W{x) < c 2 
5 / 3 W(x)) =w e (0, oo), where 



w 



poo 

A ds (2s + s 4 )" 
Jo 



(4.2) 



(4.3) 



Proof: If x' = 2ir — x, then cos y = — cos § , and sin = sin | . Therefore, 
F-j-(27r — x, 27T — y) = F±(x, y), and a change of variables shows that (PW)(x) = 
for all x. We will soon prove that the function 

r-2-K o 



f(x) = uj{x) *W(x) = u>(x)3 



dy 



o y/F + (x,y) 



(4.4) 



is continuous, with /(0) = wq > 0. This implies that / has a minimum and 
maximum on [0, 2n], and since f(x) > 0, the minimum is non-zero. This will 
directly imply that W is continuous and satisfies the bounds in (14.21) . Therefore, to 
complete the proof of the Lemma, we only need to study /. 

Suppose x £ (0, 2tv). Then the bound (13.181 ) allows using the dominated con- 
vergence theorem to prove that lining f(x + h) = f(x), which proves that / is 
continuous at x. Thus we only need to prove that / is continuous at x = and 
x = 2-7T and, as /(27r — x) = f(x) for all x, it suffices to study the limit x \ 0. As- 
sume thus < x < 7r/4. ThenforallO <y< §7r, we have cos(x/2)+cos(y/2) > 
cos(7r/8) — cos(-7r/4) > 0, and thus also F + (x, y) > C > 0. Therefore, 



lim fix) 



lim 



dy 



(4.5) 



where s x = sin(x/2), c x = cos(x/2), etc. In the final integral, let us denote 
e = s x , and change variables to s = e -1 / 3 
equal to 



sin |. This shows that the integral is 



i i 

e"3 2"2 



ds- 



4eT 



VI - e 2 /3 s 2 



Ae 1+ h + 



e 



2/3. ,2 



y/l - e 2 / 3 s 2 + Vl^e 1 



(4.6) 



Here e 2 / 3 can be canceled between the two factors. Now x — > + implies e — > + , 
and the limit can also be taken directly from the integrand, as a straightforward 
application of the dominated convergence theorem will show. Therefore, we can 
conclude (with a final change of variables to s/2) that 

4 



lim fix) 



ds 



4s + is 4 



Wq. 



(4.7) 
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Figure 1: Plot of the positive part of F_ (x, y). 



Clearly, wq is strictly positive and finite, and thus by defining /(27m) = wq, f 
becomes a function which is continuous everywhere. This completes the proof of 
the Lemma. □ 

We then require information about structure of singularities of the kernels K\ 
and K2 defined in (12.6b - K2{x, z) is bounded apart from the point singularities at 
(x, z) = (0, 2ir) and (2ir, 0), and estimate (13.181 ) will suffice to control its behavior. 
In contrast, K\(x, y) has two line singularities of strength |, which coalesce at the 
corners (x, y) = (0, 2ir) and (2tt, 0) forming a point singularity of strength 1. To 
control these singularities, we will resort to the estimates given in the following 
Lemma. For the sake of illustration, we have plotted the positive part of F_ in Fig. 
ffl 

Lemma 4.2 Let x € (0, 2tt) be given. Then there are y±,y2 such that < y\ < 
2ir — x < ?/2 < 2ir, F_(x, y) < Ofor y\ < y < y2, and 

F_(x,y) > C{yi -y)sin|, for < y < y u (4.8) 
F_ (x, y) > C(y - y 2 ) sin | , for y 2 < y < 2vr, (4.9) 

with a constant C > independent ofx,y. 

Proof: As F_(27r — x, 2ir — y) = F_(x, y), it suffices to prove the Lemma for 
< x < 7T. For notational simplicity, let c = cos | and s = sin |. Then < s < 1 
and c = Vl - s 2 G [0, 1). 
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Figure 2: Plot of d 2 F_ {x, y) (solid line) and of F_ (x, y) (dashed line) for x = 2.5. 



Since F_ is continuous and F-(x,0) = (1 + c 2 ) 2 > 0, F_(x,27r) = (1 - 
c 2 ) 2 > 0, and F-(x, 2tt-x) = -s 2 < 0, we can find < y\ < 2-ir-x < y 2 < 2ir 
such that F-(x,yi) = and F^(x,y) < for y\ < y < y 2 . Assume then that 
< y < yi, when < x + y < 2tt. As F_(x, yi) = 0, we have 



and to complete the proof of (14.8I ). it will be sufficient to show that ^^-(x, y) < 
—Cs for all < y < y\. Similarly, to prove ( 14.91 ), it suffices to show that 
d 2 F_(x, y) > Cs for all y 2 < y < 2vr. 
Let us thus consider the function 



We claim that there are y± such that F 2 (y) is strictly decreasing for < y < y'_ 
and for y' + < y < 2tt, and strictly increasing for y'_ < y < y' + . As ^(0) = 
—2s < and F 2 (2ir) = 2s > 0, then there is a unique y[ such that F 2 (y[) = 0. 
Then y'_ < y[ < y' + , F 2 {y) < for y < y[, and F 2 (y) > for y > y[. Since 
F_ is then strictly decreasing up to y = y[ and after that strictly increasing, we 
can conclude that y\ < y[ < y 2 , and that F_(x,y) > for y < y\ and for 
y > y 2 . Therefore, to complete the proof of the Lemma, we only need to find 
< C < 2 such that d 2 F^(x,y 1 ) < -Cs, and d 2 F_(x,y 2 ) > Cs. To make the 
above argument more transparent, we have plotted a sample F 2 in Fig. [2] 

Let us first consider estimating d 2 F_ (x,y 2 ). Since y 2 > 2tt — x, now cos ^ < 
— c, and thus 




y 



(4.10) 



F 2 {y) = d 2 F.{x,y) 



sin — I c + cos - — 2s cos - . 



(4.11) 




(4.12) 



Therefore, denoting t 2 = sin 




(4.13) 
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Thus if t\ < \, F 2 (y 2 ) > V2s. But also when t\ > \, F 2 (y 2 ) > 2 1 ~ 3 / 4 ^/s > 
2 1//4 s. Therefore, always d 2 F_(x,y 2 ) > 2 1//4 s. 

We estimate d 2 F- (x, y\) next. Since y\ < 2tt — x, we have cos % > — c, and 



cosy + c = 2^5 sin y. (4.14) 

Let e > be sufficiently small. If y± < tt — 2s, then cos ^- > sin e > 0, and thus 
F 2 (yi) < —2s sine. If \yi — ir\ < 2e, then | cos ^-| < sine, and sin ^- > cose. 
Thus, by (l4~T4l . 



3 



F 2 (yi) < —2s (cose) 2 — sine 



(4.15) 



Therefore, choosing e sufficiently small (for instance, s = \) we have again ob- 
tained a bound of the required type. 

We have thus proved the result for yi < vr + 2e. Assume then y\ > 7r + 2e, 
and let ti = sin ^. Then cos ^ = —\J\ — t\, and (14.141) implies that 



2vtoi = c-Vl-tt = 1 < -^==- (4-16) 



Therefore, t^ 2 > 4-^/1 — tf^/s, and thus 



F 2 (yi) = -2bUI + 2syi - tf < -2sWl - if (4 - 1) < -6ssine. (4.17) 

This proves that there is a pure constant C > such that d 2 F^ (x, y\) < —Cs. 

We still need to prove the monotonicity property of F 2 mentioned earlier. Let 
u = cos |, when u goes from 1 to —1 strictly monotonicly, as y goes from to 2-k. 
Also 



F 2 (y) = -2su - (c + u) a/I - ^ 2 = s(tt)- (4.18) 

Then 



V \ o /i 2 , n ( c + ' u ) n, \ c + 3u- 2u 3 

,(n) = -2s- V / l-^ + - 7 _, g («) = (j _ ^ . 



(4.19) 



The polynomial c + 3u — 2u 3 has a local minimum at u = —2~ 1 / 2 and a local 
maximum at u = 2 -1 / 2 . Its values at u = —1, 0, 1 are c — 1, c, c + 1, respectively. 
Thus there is — 1 < u$ < such that g"(uo) = 0, g"(u) < for u < uq and 
g"{u) > for u > uq. Since c < 1, g'(u) first decreases strictly from +oo to 
g'(uo) and then increases strictly to +oo again. Since g'(0) < 0, also g'(uo) < 
and thus there are u± such that g is strictly increasing for — 1 < u < u_ and for 
u + < u < I and strictly decreasing for u_ < u < u + . This implies the stated 
monotonicity property of F 2 and completes the proof of the Lemma. □ 
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We then prove two intermediate compactness results which will become use- 
ful in the proof of Proposition 12.41 Instead of Sobolev-space techniques (such as 
proving that the operators improve a Sobolev index), we will rely on direct norm 
estimates which are quite straightforward in the present case. 

Proposition 4.3 Let ip : / — > C be given, and assume that there are C,p > such 
that 

\i/)(x)\ <c(sin§Y (4.20) 



2 

for all x £ I. Then the function 

K(x,y) = t/j(x)*K 2 (x,y)^(y) (4.21) 
defines a compact, self-adjoint integral operator on L 2 (I). 

Proof: Since (sinf) p ~ 1/2 e L 2 {I), the estimate (l3~T8T) proves that K is Hilbert- 
Schmidt, and thus also compact. As K(x, y) is symmetric, the operator is self- 
adjoint. □ 



Proposition 4.4 Let if) : I — ► C satisfy the assumptions of Proposition \473\ Then 
the function 

K(x,y) = il>(x)*K 1 (x,y)xlj(y) (4.22) 
defines a compact, self-adjoint integral operator on L 2 (I). 

Proof: Let p, C > be constants for which (14.201) holds. As the bound is a de- 
creasing function of p, it is sufficient to prove the Proposition assuming < p < I. 
Suppose < e < 2tt is arbitrary, and let T £ denote the integral operator defined by 

X £ = {(x,y)e [e,2n-e] 2 \F_(x,y) >Cesin|}, (4.23) 
T e (x,y) = t((x,y) G X £ )K(x,y). (4.24) 

Then \T e (x, y)\ < c(esin I)" 1 / 2 for some constant c, and since the kernel T £ (x, y) 
is obviously symmetric, we can conclude that T £ is a self-adjoint Hilbert-Schmidt 
operator on L 2 (I). As K(x,y) is also symmetric, we can apply Proposition IA.1I 
to the integral operator K — T £ . We will choose cj){x) = sin(a;/2) p , a = 1 — 
i and A(x,y) = sm(x/2)~P(K(x,y) - T £ (x, y)) sin(y/2)-P. By (l4~20l) . then 
|^4(a^,y)| < C 2 l((x, y) X £ )K\(x,y). Therefore, it is enough to inspect the 
integral 

J(x) = jd y (sm|) P+ ' (sm|) P ~ 5 l((x, y) X^K^x, y). (4.25) 
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We claim that there is c > 0, such that J(x) < ce p for all x € I. Then by 
Proposition IA.1[ we have \\K — T e \\ < C 2 ce p . Since this also holds for e > tt, 
when T £ = 0, we find that K is itself a self-adjoint bounded operator. However, 
then also T £ — > K in norm, and since each T e is a compact operator, and the space 
of compact operators is closed in the operator norm, we conclude that K is also 
compact, proving the results mentioned in the theorem. 

Thus we only need to show that J(x) < ce p for all x € I, assuming < p < I. 
For any iGlwe obtain from Lemma l4~2l the following rough estimate, where the 
integration region is estimated trivially, 

J(x) < C- 1 ' 2 (sin |) P ( jT dy h(y) + J J dy f 2 (y)) , (4.26) 

where 

f^y) = Un V -Y~~ 2 — = and f 2 (y) = (sin y -Y~~ 2 — =. (4.27) 
v 2/ Vyi - y V 2/ ^y - y 2 

Since 

f'Vi /"l 

/ dy^-5( yi -y)-2 =y P / dtt p -5(l-t)-5, (4.28) 
io Jo 

both of the remaining integrals are uniformly bounded in x, independently of the 
actual values of y\ and y 2 . The rough estimate proves that sup^ J(x) is uniformly 
bounded for all e, but it also proves that if x S [0, e) or if x £ (2-7T — e, 2ir), then 
there is a constant d such that J(x) < c'e p . 

Let us then consider the remaining case when e is small (say e < 1) and x G 
[e, 27r — e]. Then sin | > sin |, and thus Lemma l4T2l shows that every (x, y) E I 2 
for which e <y < y\ (x) — e or y2(^) + e < y < 27r — e, belongs to X e . Therefore, 
such y do not contribute to J(x), and we can estimate 

J(x)<C-V 2 Un -) (/ dyf 1 (y)+ dy h(y) 

v /7 V J0 Jmax(0,i/i-e) 
/•min(27r,j/2+e) /"27r . 

+ / dy/ 2 (y)+ / dy/ 2 (y)J. (4.29) 

^ ?/2 max(27r— e,y2) 

Each of the four integrals can be estimated similarly to (14.281) . which shows that 
they are bounded by c"e p for some constant c". Therefore, we have proven that also 
in this case J(x) < ce p for some c. This completes the proof of the Proposition. □ 

Proof of Proposition \2.4\ By Propositions 14.31 and 14.41 the operators ioK\io and 
luK 2 u) are compact, and thus also bounded. By Lemma 14- 1 1 W is a bounded 
function, and thus L is a sum of a bounded multiplication operator and a compact 
integral operator. Thus L is bounded, and then the argument in the beginning of 
the section, based on Proposition 13. 1[ implies that it is positive. On the other hand, 
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Lemma |4j] also implies that < V(xy 1 / 2 = oj(x)W{x)- 1/2 < c 2 sin(x/2) 1 / 6 , 
and we can apply Propositions 14.31 and 14.41 also to the definition of B, equation 
(12.91) . This proves that B is a compact, self-adjoint operator on L 2 (I). B and L 
also commute with P, by the symmetry properties of F± stated in the beginning of 
the proof of Lemma |4~T1 

Thus we only need to prove the last claim in the Proposition. Applying the def- 
initions, we find L = W l l 2 (l - B)W 1 / 2 . For any ib G L 2 , for which W~ l / 2 ib G 
L 2 (I), we can find a sequence f n of Lipschitz continuous functions, such that 
f n — > W~ l / 2 ib in L 2 . Thus by the boundedness of L and Proposition 13. 1[ then 

<</>, (1 - B)if}) = lim(/ n , Lf n ) = lim / dxdz JL = 

x |flr n (ar) + g n {h(x, zj) - g n (z) - g n {x - z + /i(x, z))\ 2 , (4.30) 

where g n = w/ n — > F _1//2 ^. The function defined by the integral on the right 
hand side is L 2 -continuous in f n . To see this, let us inspect the difference of two 
such integrals, which can be bounded by a sum of finitely many terms of the type 
Jdxdz{4F + )-^ 2 \G(X)\ 2 , where X denotes any one of the functions x, z, h(x, z), 
or x — z + h(x, z), and G is in L 2 . The first two choices of X lead to integrals 
which clearly can be bounded by J dxV (x)\G(x)\ 2 . However, so do the last two 
choices, as can be seen by employing the symmetry h(z, x) = x — z + h(x, z) 
and Lemma 1331 Using then the fact that any relevant G is of the form G = ujF, 
F G L 2 , we have here J dxV{x)\G{x)\ 2 = J dxW{x)\F{x)\ 2 < HWHooH-FH 2 . 
This suffices to prove the continuity, and thus for the above class of if), 

lib, (1 - B)ib) = / dxdz — - 1 

JP 2^F^z) 

x \g{x) + g(h(x, z)) - g(z) - g(x- z + h{x, z))\ 2 , (4.31) 

with g = V~ l l 2 ib. Then the previous argument can also be applied to show that 
the right hand side is L 2 -continuous in ip, which proves that (14.311 ) holds for all 
ib G L 2 . Therefore, 1 — B > 0, and Bib = ifi if and only if the integral on the 
right hand side of (14.311 ) vanishes for g = V~ x l 2 i\). Since then the integrand must 
be zero almost everywhere, this is possible if and only if the periodic extension of 
g is a collisional invariant. □ 



5 Collisional invariants (proof of Theorem 12^2]) 

It is clear that every ib(x) = c\ + C2U>(x) is a locally integrable collisional invari- 
ant. Thus to prove the Theorem, it will be enough to consider any ip, which is 
a locally integrable collisional invariant, and to show that it is almost everywhere 
equal to a function of the above form. Let us assume ib is such a function. Then, 

as Q(x, h(x, z), z) = for all x, z, we have for almost every i,z£M 

ib(x) + ip(h(x, z)) — ib(z) — ib{x — z + h(x, z)) = 0. (5.1) 
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In addition, since h(2ir — x,2tt — z) = —h(x, z) for x, z G /, then also Pip sat- 
isfies (15.11 ) almost everywhere. Therefore, both of (ip ± Pip)/ '2 have this property, 
and thus it is sufficient to prove the result assuming that ip is either symmetric of 
antisymmetric under P. 

Let us begin by showing that then there is / : R — ► C which is periodic and 
twice continuously differentiable apart possibly from points in 2ttZ, and for which 
ip = f almost everywhere. We will do this by integrating (15.11 ) over x. However, 
the integration region has to be chosen with some care, in order to guarantee that 
the result is finite. With a certain < Eq < f to be fixed later, we consider an 
arbitrary < e < Eq. We define for all e < z < tt + e 

1 r 27T 

f h£ (z) = - dx[iP(x)+iP(h(x,z))-iP(x-z + h(x,z))]. (5.2) 

£ J2-n-e 

A comparison with (15.11 ) reveals that then fi )£ (z) = ip(z) almost everywhere. On 
the whole integration region x > z and thus 

z — X x — z 
h(x,z) = — \-$>(x,z) and x — z + h(x, z) = — \-<S?(x,z) (5.3) 

with 

/ x — z £ + 
<P(x, z) =2 arcsin I tan — - — cos — - — J . (5.4) 

We claim that there are Eq,C > such that for all above 

2d x $(x,z) < -1-C. (5.5) 

Together with (15.31 ) this implies that the last two mappings in the arguments of tp 
in (I5.2I ) are strictly decreasing in x. In particular, when x — > 2ir, we have also 
h(x, z) \ z — 2-k and x — z + h(x, z) \ 0. Thus a change of variables and 
denoting Co = h(2ir — e, z) yields 

1 r 2n 1 f£o 9 

f hE (z) = - dxiP(x) + - dC^(C) 



£J2k-s £Jz-2tt l-2d x ${x 1 ((,z),z) 

I /"STT-e-z+Co 



e J -l-2d x <S>(x 2 ((,z),z) 

Since both of the factors multiplying ip(C) WQ continuous in z and uniformly 
bounded, we can conclude using the dominated convergence theorem that /i e is 
continuous. 

As mentioned before, ip{z) = fi >E (z) for almost every e < z < it + e. Since 
Pip = aip, with a G {±1}, we can define f £ (z) = fi, £ {z) for e < z < tt + e, 
and f £ {z) = afi !£ (27r — z) for tt — e < z < 2ir — e. In the common domain 
near z = tt both functions have to be equal everywhere, as they are continuous 
and coincide with ip almost everywhere. In particular, f £ (z) is also everywhere 
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continuous and equal to ip almost everywhere. Since e was arbitrary, we can then 
extend the definition to cover the whole of (0, 2tt), by choosing f(z) = f e (z) for 
any e < z, 2ir — z. Again, by continuity, any two functions f £ and f £ i must agree 
on the intersection of their domains of definition, so / is a continuous function on 
(0, 2tt), which we extend periodically to K. Then ip(z) = f(z) a.e. z£l. 
Using the continuity of h, this implies that 

f(x) + f(h(x, z)) - f(z) -f(x-z + h(x, z)) = 0, (5.7) 

for all x, z £ I for which all arguments are non-zero, i.e., whenever x ^ 0, z ^ 
and i/z. In particular, then (15.61 ) holds for all e < z < tt + e after both ifi and 
fl >£ are replaced by /. However, then the right hand side of (15.61) is continuously 
differentiable, and we can conclude that / is continuously differentiable on (0, 27r). 
This argument can then be iterated once more to conclude that / must be twice 
continuously differentiable on (0, 27r) (this way even smoothness could be proved, 
but we will not need this property here). 

We next prove that we can choose /(0) so that / is continuous and f'(x) has a 
limit for both x \ and for x / 2ir. Since h(x, 2ir — x) = tt — x for all x G /, 
we have for all x G (0, tt), 

f(x) - /(2vr -x) + /(tt - x) - /(tt + x) = 0. (5.8) 

If / is antisymmetric, f(ir) = and (15.81) implies that for all x G (0, tt), f(x) = 
/(•7T + x). Therefore, in this case / is continuously differentiable at x = 0, after 
we define f{2irn) = 0, n G Z. 

Assume then that / is symmetric which implies f'(z) = —f'(2ix — z). Let us 
consider values < z < x < 2ir, when (15.31 ) holds. Differentiating (15.71) with 
respect to x and z yields 

f(x) + d x hf'(h) - (1 + d x h)f(x -z + h)=0, (5.9) 
-/'(*) + d z hf'(h) - (-1 + d z h)f'(x -z + h) = 0. (5.10) 

We multiply the second equality by (l+d x h), and then use the first one to eliminate 
f'(x — z + h). This proves that 

(i - d z h)f(x) - (i + d x h)f( z ) + (d x h + d z h)f(h) = o. (5.ii) 

We divide the equality by 1 — d z h and consider taking the limit x — ► 2ir for a fixed 
z. Then h — > 2 — 27r, and the partial derivatives converge as (see (15.171 ) to obtain 
explicit formulae from which these can be checked) 

d x h(x, z) -» -(1 + t 2 ) and 9 z /t(x, z) -> 1, (5.12) 

where t = tan 2?r ~ z . Since 1 — d z h — * 0, we need to compute the limit more 
carefully. Let us fix for definiteness, z = tt, when a straightforward computation 
shows that d z d x h{x, z) — > 1, and thus by LHospital's rule, 

Hm f'(z)-f'(h(x,z)) = ^ f"(h(x,z))d x h(x,z) = _ 2f „ {n) (5 13) 

x/"27r l — d z h(x,z) x d x d z h(x,z) 
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We can then use this result in (15.1 lb to prove that the limit of f'(x) exists when 
x /* 2ir, and thus by symmetry the same is true about the limit when x \ 0. 
Since this implies that /' is bounded on [0, 2tt}, we can also conclude that the limit 
c = linx^o / 0*0 exists, and we can make / continuous by defining /(27m) = c 
for all n G Z. 

We can thus assume that / is continuous and periodic on R, continuously dif- 
ferentiable on (0, 2tt), and that a = linx^o f'( x ) an d b = ]iva x fix) exist. 
We also have b = —a, if / is symmetric, and b = a, if / is antisymmetric. Let us 
now consider any < z < it, and x = 2tt — e for < e < 2tt — z. As proven 
earlier, in the limit e \ 0, x — z + h\ 0, and we get from (15.71) that when e — > 0, 



/(2tt + h)- f(z) f( x + h-z)- /(0) + /(2tt) - f(x) 



2tt — x 2ir — x 

-» -(1 + d x h)a + b. (5.14) 

The left hand side converges to —f'(z)d x h, and, since by (15.121 ) d x h(2TT,z) = 
— (1 + t 2 ) < 0, we have proven that 

f'(z) = a- T ^(a-b), (5.15) 

where = cos 2 = | (l — cos §) = | — 0/(2). Therefore, for < z < 7r 
we need to have 

f'(z) = ^ + (a-b)u'(z). (5.16) 

If / is antisymmetric, then f'(z) = a for < z < ir. Since then also /'(27T — 
z) = a, we must have f(z) = c + az for < z < 2-k. However, as also /(0) = 
= /(27r), we need to have c = = a, and thus the only antisymmetric solution 
is the trivial solution / = 0. If / is symmetric, b = —a, and /'(27r — z) = —f'(z). 
Thus then f'(z) = 2auj'{z), for < z < 2ir, and there is c 6 C such that 
f(z) = c + 2auj(z) for all z£K. Therefore, / is in both cases a trivial collisional 
invariant, and since ip = f almost everywhere, we have arrived at the conclusion 
made in the Theorem. 

We still need to prove (I5.5I ). Using the shorthand notations t = tan c = 
cos ^^jp, we can write 



< c +\=* , „ y ) (5 .,7, 



vi - c 2 t 2 V 1 + * 2 t\/r^ 2 + vi - c 2 t 2 

For e, x, z as above, i.e., for0<e< : |,e<2;<7r + e, and 2n — e < x < 2ir, 

7T E X — Z 7T £ 7T Z — E X + Z 3n £ 

< < and — I < < h-. (5.18) 

4 2-4-2 4 2 4-4-4 4 V 7 
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Thus t > tan (f - §) > 0, and with d = cos ^ 

1< -d < c < cos ( - + ) = -sin^-^ < 0. (5.19) 

" ~ 1 2 4 / 4 ~ 

If also z < ^ — e, then x — z > ^ and thus t > tan(37r/8) > 2. In this case, we 
can estimate the first term using c < 0, which yields 

t 2 - 1 

2d x $(x,z) + l<--—- = -(t-l)<-l. (5.20) 
t ~r J. 

Otherwise, z - e > § - 2e > § . If t > 1, we find 

2<9 1 .$(x, z) + 1 < c < - sin < - sin — < 0. (5.21) 

4 16 

On the other hand, if t < 1, then by 



1 — t 2 X — Z /7T 



cos < cos e = sine, (5.22) 



we have 



1 + t 2 2 ~ \2 



2d x $(x,z) + l< , o (c + 



VI - C 2 t 2 V 1+t 2 ^/l - ( C ')2 



1 + t 2 / . 7T sin e 



^ /i 272 " Sin Tfi + • ( 5 - 23 ) 

VI - c 2 i 2 V 16 sin 

Since the term in the parenthesis approaches — sin < 0, when e — > 0, there is 
eo > such that for all < e < eo the right hand side is less than — | sin Thus 
we can conclude that for all such s Equation (15.51 ) holds at least with C = \ sin ^ . 
This completes the proof of the Theorem. 



6 Resolvent expansion (proof of Theorem 12.51) 



Let us begin with the following corollary of the results proven in the previous 
sections. 

Corollary 6.1 The eigenspace of B with eigenvalue 1 is two-dimensional, and it 
is spanned by V(x) 1 ^ 2 and uj(x)V {x) l l 2 . Every ip S L 2 (I), such that Pxfj = —if), 
is orthogonal to this eigenspace. 

Proof: Suppose tp belongs to the above eigenspace of B. By Proposition 12.41 then 

is a collisional invariant. Thus by Theorem 12.21 there 
are c\ , C2 G C such that cuW^ 1 ^ 2 ^ = c\uj + C2- By Lemma l4~Tl both of the vectors 
W l l 2 = ujV 1 / 2 and uo~ l W l l 2 = V 1 / 2 belong to L 2 and are symmetric under P. 
This proves the results stated in the Corollary. □ 
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To estimate error terms, we will rely on the following estimates: 
Lemma 6.2 There is C > such that for any < A < 1, 

^dx T - fsin;r*<CA-*. (6.1) 



W(x) + A V 2 

Proof: Let < A < 1 be arbitrary. By symmetry, the integrals over [0, n] and 
[tt, 2tt] are equal. On the other hand, by Lemma |4~T1 we have 



Jo 



'dx "& (0n^<-±- fdx jf . (6.2) 



The integral over x G [^7r, 7r] is clearly bounded uniformly in A. To estimate the 
integral over [0, \tt], we change the integration variable to s = X~ 3 ^ 5 s x , which 
shows that 

/•tt/2 s 1 / 2 3 3 f°° s l/2 t 

/ dx — I <2^Ai+to- 1 / ds <cA-io, (6.3) 

Vo 4 7 +A Vo s 5 / 3 + l 

since the integral over s is finite. Thus (16.11 ) holds for some finite C. □ 
Lemma 6.3 For < A < 1, let 

^ = B w^x J - (6 " 4) 

For any < e < ^ there is a constant c £ > such that for all x £ I amJ A, 

Iv?a(x)| < c £ A-^- £ (v(s) sin |) ~ 5 . (6.5) 

Proof: Let < e < be arbitrary. Applying the definitions of B and W, as well 
as the bound \u>'\ < i, we find 



|v>a(*)| < l -V{xY 1 /* j\y {\K l {x,y)\+2\K 2 (x,y)\) w ^ ] +x . 
By Lemma [6721 and Eq. (13- 18b . the second term in the sum satisfies 



(6-6) 



" dy2|K2( ^ y)l ll^^( Sm f)" lA ~"' (6 " 7) 
and thus leads to a bound of the desired form. 
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To the first term we apply Lemma l4~2l which shows that 

< C 



fc *|jr l( «,,)| 



dy 



W(y) + A 
1 



sm ■ 



+ 



2tt 



dy- 



vW^ + A -W) y / y-V2(x) W{y) + A 

(6.8) 



Changing the integration variable to y' = 2ir — y in the second integral reveals 
that it is equal to the first integral, if y\{x) is replaced by 2-k — y^ (x). Thus it is 
sufficient to inspect the first integral. We estimate it using Lemma l4~T1 and Holder's 
inequality for p' = 34 _^ 0g < 2, q' = 3 _ 6 10£ > 2. This shows that for any y\ G /, 



dy 



3+lOe ^ 

«>(j/) 



< 



l + Cl 



.yi 



dy(yi - y) 2 



1 













w(y) 5 / 3 + A 



(6.9) 



The first factor is an e-dependent, finite constant, and the second factor can be esti- 
mated as in the proof of Lemma loT2l first we use symmetry to reduce the estimate 
to [0, 7r], then u)(y) > on [7r/2,7r] to bound the integral on this interval by a 
constant, and finally on the interval [0,tt/2] we change the integration variable to 
s = A~ 3 / 5 sin |. This shows that there is a constant c such that 



2tt 



dy 



w(y) 5 / 3 + A 



< ca! ^'(f- 1 ) 



ds 



S 5 /3 + 1 



<c>\-h'+l, 



(6.10) 



where the integral over s is finite, since q' > 2 > |. Therefore, as — | + ^ 
— ^ — e, (16.91 ) implies that (16.51 ) holds also for the first term in the sum. 

This has the following immediate corollary: 



□ 



Corollary 6.4 For any < e < | there is a constant c £ > such that ||<^a|| 2 < 
c e A~5~ £ for all < A < 1. In addition, P(f \ = —<p\- 

Proof: Let e' = e/2. Then by Lemma [531 

-2tt . r 2n 



\<P\\ 



dx \ ip\(x)\ < c e i\ 5 



-2e' 







dx ( V(x) sin — 



(6-11) 



By Lemmal4Tl there is C such that V(x) sin(x/2) > C(sin(x/2)) 2 / 3 , and thus the 
remaining integral over x is finite. This implies that there is c £ such that ||(^a|| 2 < 
c £ \~5~ £ . Since both B and W commute with P, and u>' is antisymmetric, it 
follows that ipx is antisymmetric. □ 
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Armed with the above results, we can prove the main theorem. We claim that 

for any < e < f , 

, 1 



R(X) 



x + w 



uo')+0{\- 



(6.12) 



This implies that only the first term of a resolvent expansion needs to be considered 
for the limit (12.131) . However, then 



LO 



X + W 



-lo 



2tt 



dx 



a! 



2 ^ 

A + W(x) C ° S 2 



(6.13) 



and we can use the symmetry of the integrand to reduce the integration region to 
[0, 7r], while gaining a factor of 2. We then change variables to s = A~ 3 / 5 sin |. By 
Lemma |4~T1 the remaining integrand is dominated by 1+e * s5 / 3 which is integrable 
on (0, oo). Thus 



lim (a/, — 

A ^o+ \ X + W 



LO 



CO 



Jo 1 + U 



5 • 
+ WoSS 



(6.14) 



Clearly, < cq < oo, as claimed in the Theorem. 

To estimate the "error term" in (16.121) assume < A < 1 is given. Let A = 
W^BWz, when A is a bounded operator and L = W — A. We use the resolvent 
expansion of L up to the second order, 



1 



1 



+ 



1 A 1 1 A 1 A 1 

— + —A ^A- 



X + L X + W X + W X + W X + W X + L X + W 



(6.15) 



where, since W > 0, (A + W) 1 is a bounded, positive operator. Therefore, 
denoting (f> x = A^w', 



(6.16) 



Using Proposition 12.41 Corollary 16-11 and Lemmas I6.2146.4i we can now prove 
(16- 12b - To estimate the first correction we use Lemmas 16. 21 and [631 



LO 



< 



1 :A'„> 



X + W X + W 



W2 

w + x 



^ dx „^ X \ . V(x)$\<px(x)\ < c e CX~ l s- 



o 



W(x) + X 



(6.17) 



This proves that the first correction is of the claimed order, and we only need to 
inspect the final term in (I6.16I ). 
Firstly, 4>\ = W^ifx, and thus 

1 . \ / 1 



X + L 



1-B + XW- 1 



(6.18) 
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Here, by Lemma |4~T1 we have the operator inequalities 

1 -B + XW~ l > 1 -B + — > — > 0. (6.19) 

where ^ is proportional to the unit operator, and thus commutes with B. There- 
fore, 

(l-B + XW- 1 )- 1 < (1-B + — y 1 (6.20) 

implying 

£ (*»• 3TTT* A ) - (**• i + v« <6 ' 21) 

By Proposition 12 .4 1 is compact, and its spectrum consists of isolated eigenvalues 
(apart from zero). Since also B < 1, we can order the eigenvalues so that 1 = 

Ai > A2 > In particular, then 5 = 1 — A2 > 0. Since tp\ is antisymmetric, it is 

orthogonal to the eigenspace of B with eigenvalue 1 by Corollary 16.11 Therefore, 
using the spectral decomposition of B, we find that 

(^■ i +v !»-fl yA )*?+W w "'-^~ t ~'' <U2) 

where we used the estimate in Corollary 16.41 Then we can conclude from (16.211 ) 
that (16.121 ) holds. This completes the proof of Theorem 1231 



A Integral operators 

Given a positive measure n on X, any function K : X x X —> C, which is 
measurable in jj, X \i, can be used to define an operator T in L 2 ((i) by the formula 



(Tf)(x)= / n(dy)K(x,y)f(y). (A.l) 
Jx 

More precisely, we define T as a possibly unbounded operator with the domain 

D(T) = { f G L 2 (fi) [ n(dx) ( [ fi(dy)\K(x,y)\\f(y)\) < 00} . (A.2) 



K is then called the integral kernel of the integral operator T. 

We need here only the following convenient estimate for an operator norm of 
such integral operators. 

Proposition A.l Let fi be a positive measure on X, and assume that A : X x 
X — ► C is measurable with respect to p, x p, and satisfies A(x, y)* = A(y, x) for 
almost every (x, y) £ X 2 . Consider any measurable (j> : X — > C, let B(x, y) = 
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(f>(x)*A(x, y)4>(y), and let T denote the corresponding integral operator. If there 
exists a£i such that 



<7 a = esssup I \4>{x)\*- a j p(&y)\A(x iy ) \ \4>{y)\ a j < oo, (A.3) 

then T is a bounded, self-adjoint operator on L 2 (fi), and \\T\\ < C a . 

Proof: Let a G M. be given, and let us denote a' = 2 — a. Then |^(a0| 2 = 
|(/>(rE)| a |(/)(x)| Q . (We apply here the usual convention used in connection with 
positive measures, that • oo = 0. Let also 0° = 1). Thus for any / £ I 2 , we have 
by the Schwarz inequality and Fubini's theorem 

(jjl x fi x fj,)(d(x,y,z))\f(x)\ \B(z,x)\ \B(z,y)\ \f(y)\ 

(fix fix (i)(d(x,y,z)) 



x \f(x)\ [\A(z,x)\ \A(z,y)\ \<l>(x)\ a \^{z)\ a+a 
x \f(y)\ (\A(z,x)\ \A(z,y)\ |^(»)| a '|^)| 0,+0, '|^)r 
< / Md*)|/(*)| 2 (V(*)r' / lx{dz)\A{x,z)\\4>{z)\ a 



x \<t>(zT J Kdy)\A(z,y)\\4>(y)\ a 

< C 2 a \\f\\ 2 < oo, (A.4) 

where we have used the symmetry of A. Therefore, D(T) = L 2 , and since the left 
hand side of (IA.4b is an upper bound for ||T/|| 2 , we have proven that ||T|| < C a . 
As B{x, y)* = B(y, x) almost everywhere, T is then also self-adjoint. □ 
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